HIGHER ORDER ELLIPTIC AND PARABOLIC SYSTEMS 
WITH VARIABLY PARTIALLY BMO COEFFICIENTS IN 
REGULAR AND IRREGULAR DOMAINS 



HONGJIE DONG AND DOYOON KIM 

Abstract. The solvability in Sobolev spaces is proved for divergence 
form complex-valued higher order parabolic systems in the whole space, 
on a half space, and on a Reifenberg flat domain. The leading coeffi- 
cients are assumed to be merely measurable in one spacial direction and 
have small mean oscillations in the orthogonal directions on each small 
cylinder. The directions in which the coefficients are only measurable 
vary depending on each cylinder. The corresponding elliptic problem is 
also considered. 
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1. Introduction 

We study the solvability in Sobolev spaces for parabolic operators in di- 
vergence form of order 2m 

Ut + {-irXu:=Ut + {-ir Yl D^ia^pD^u), (1.1) 

a|<m,|/3|<m 

where a and /3 are multi-indices, a^^ = [a^^(x)]"j^;^ are n x n complex 
matrix- valued functions, and u is a complex vector- valued function. As 
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usual, for a = (ai, . . . , a^), we write 

All the coefficients are assumed to be bounded and measurable, and L is 
uniformly elliptic; cf. ()2.ip . In the case that all the coefficients and functions 
involved are independent of the time variable, we also deal with elliptic 
operators in divergence form of order 2m as in (jl.ip without the Ut term. 

Our first focus in this paper is to find minimal regularity assumptions on 
coefficients for the Lp-solvability of higher order elliptic and parabolic sys- 
tems. In other words, we prove that there exist unique solutions in Sobolev 
spaces (e.g., Ty^(Q), 1 < p < oo in the elliptic case) to higher order elliptic 
and parabolic systems with coefficients having less regularity assumptions 
than those available in the literature. We call the class of coefficients in the 
paper variably partially BMO coefficients, the key feature of which is that, 
on each cylinder (or ball in the elliptic case), the coefficients are allowed 
to have no regularity assumptions (i.e., merely measurable) in one spatial 
direction. If we name the spatial direction the measurable direction, variably 
partially VMO coefficients mean that, for each small cylinder, there exists 
a measurable direction, which may depend only on the cylinder, such that 
coefficients are allowed to be merely measurable in that direction and have 
small mean oscillations in the orthogonal directions. See Section [2] for a 
precise formulation of the assumption. This class of coefficients was first 
introduced by N. V. Krylov about two years ago in [29] in the context of 
second order elliptic equations in non-divergence form. The same class of 
coefficients is studied in [15] for second order elliptic and parabolic systems 
in divergence form defined in the whole space. See also a more recent paper 
[6] for second order elliptic equations without lower-order terms and with a 
symmetric coefficient matrix on a bounded Reifenberg flat domain. 

In order to put into perspective the class of coefficients in this paper, let 
us review very briefly the coefficients for Lp-theory of elliptic and parabolic 
equations/systems studied in the literature. A rather rich collection of pa- 
pers about Lp-theory is devoted to the study of equations with constant 
or uniformly continuous coefficients. For instance, see [3l (U [33] for ellip- 
tic systems and [371 EOl US] for parabolic systems. There are also a large 
number of papers concerning equations with discontinuous coefficients. An 
important example of discontinuous coefficients is a family of functions with 
vanishing mean oscillations (VMO), which was firstly considered in [9| [71112] 
for second order elliptic and parabolic equations. With this class of coeffi- 
cients, Lp-estimates for higher order equations have been obtained in [10] . 
[21] , [34] , [35] and [32] . We also refer the reader to [31] for elliptic equations 
with piecewise constant coefficients, and [TT] for results of elliptic equations 
with measurable coefficients for a restricted range of m or p, and [191 E] for 
Holder estimates of systems with uniformly continuous or VMO coefficients, 
as well as [5] for a result of fourth order elliptic equations in nonsmooth 
domains. 
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Recently, in [16] we considered higher order systems in both divergence 
and non-divergence form with coefficients having locally small mean oscil- 
lations with respect to the spatial variables (and measurable in the time 
variable in the parabolic case). We call coefficients in [16] BMO coefficients, 
the class of which includes VMO coefficients as a proper subset. In this pa- 
per the coefficients also have small mean oscillations, but need no regularity 
assumptions in one spatial direction (determined locally). Thus especially 
in the elliptic case the coefficients here are strictly more general than those 
studied in [16]. For the parabolic case, once a measurable direction is de- 
termined on each cylinder, the coefficients are required to have small mean 
oscillations in the orthogonal directions, which include the time direction, 
whereas the coefficients in [16] are merely measurable in time, but have 
small mean oscillations in all the spatial directions. Thus the two classes of 
coefficients here and in [16] are not mutually inclusive. 

Our second focus in this paper is to examine the Lp-solvability of elliptic 
and parabolic systems in various domains. We first deal with systems in 
the whole space, where interior Lp-estimates are the key ingredients. Then 
we prove the solvability of systems defined on a half space, where we obtain 
boundary Lp-estimates. As usual, combining the interior and boundary Lp- 
estimates yields the solvability of systems defined on a domain as long as the 
boundary is regular, for instance, if the domain is a Lipschitz domain with a 
small Lipschitz constant. Indeed, in [16] elliptic and parabolic systems with 
BMO coefficients are studied in the whole space, on a half space, and on 
a Lipschitz domain. Here we further examine the Lp-solvability of systems 
with variably partially BMO coefficients defined on Reifenberg fiat domains 
(possibly unbounded), which are more general than Lipschitz domains. As 
is well known, at small scales the boundary of a Reifenberg flat domain 
can be approximated by hyper-planes. We remark that, for the Reifenberg 
flat domain case, on each small cylinder centered at the boundary, we as- 
sume that coefficients are merely measurable in the direction normal to the 
approximating hyper-plane. 

In the half space case, compared to previous results in \22\ [23| [14] for 
second order equations, one of the virtues of our results is that on each 
small cylinder centered on the boundary, we do not require the measurable 
direction to be exactly perpendicular to the boundary, but allow it to be 
sufficiently close to the normal direction. To this end, our proof is founded 
on a delicate cutoff argument combined with a higher order Hardy-type 
inequality; see Lemma 17.91 Indeed, due to the nature of the boundary of 
Reifenberg flat domains, studying the half space case with almost normal 
measurable directions is readily applied to obtaining the key estimates for 
systems on Reifenberg flat domains with variably partially BMO coefficients. 
For further discussions about equations/systems on Reifenberg flat domains, 
see Section [H 
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The results in this paper and in [16] imply several generalizations of the 
interior Holder estimate for higher order elliptic equations with VMO co- 
efficients proved in [2] Proposition 50. Indeed, by the Sobolev imbedding 
theorem, we obtain both interior and boundary Holder estimates for higher 
order elliptic and parabolic systems with coefficients in [IB] and with coef- 
ficients considered in this paper. 

It is worth mentioning that we impose the uniform ellipticity condition 
in this paper, while the main results in [16] are established under a slightly 
weaker Legendre-Hadamard ellipticity condition. Because the coefficients 
are measurable in one spatial variable, when getting the -L2-estimates (or 
equivalently the Garding inequality) the Fourier transform method used in 
|16j is not applicable here. This is the only place where the uniform ellipticity 
condition plays its role. 

Our proofs for the whole space and half space cases are in the spirit of 
|27| . in which the author gave a novel method of studying Lp estimates of 
second order elliptic and parabolic equations in the whole space with rough 
coefficients. Differed from the arguments in [9l [121 [lOl [35] , which are based 
on singular integrals and commutator estimates, the crucial step in [2^ is to 
establish certain interior mean oscillation estimates of solutions to equations 
with 'simple' coefficients, which are measurable in some directions and inde- 
pendent of the others. As a consequence, VMO coefficients are treated in a 
rather straightforward way. The method in [27] was later further developed 
in a series of papers |28l|22l[23l[2ll[25l[29l[IZl[ia[l3] on Lp-estimates of 
second order equations with rather general discontinuous coefficients. 

One of the main differences in our study between second order equations 
and higher order equations is the following. As is explained in [221 [23l [E] , 
when dealing with second order equations (or systems) on a half space with 
homogeneous Dirichlet or Neumann boundary conditions, we used the odd 
and even extension techniques by observing that the odd or even extension 
of the solution satisfies the equation extended in the whole space. Then we 
were able to rely on the fact that the coefficients are allowed to be merely 
measurable in one spatial direction. Thus the boundary Lp-estimates were 
derived almost immediately from interior estimates. However, the extension 
techniques do not work for higher order equations or systems, since in gen- 
eral extensions over the whole space do not satisfy the extended equations 
and, in the non-divergence case, they do not even belong to the correct so- 
lution spaces. Because of this, in this paper we use a modified version of 
a generalized Fefferman-Stein Theorem developed in [29] in order to pro- 
duce boundary mean oscillation estimates of solutions. This approach for 
boundary Lp-estimates is applicable to a wide class of equations or systems. 

Another worth noting difference comparing to second order equations is 
in the proof of interior Holder estimates; see Lemmas 14.51 and 14. 6[ To get 
around the lack of regularity of Diu, the main idea in [15] is to estimate 
instead the interior Holder norm of certain linear combination of derivatives 
of the solution which in particular contains DiU. For higher order systems. 
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the proof is more involved and we also make use of some Lp estimate of a 
one dimensional problem; see Lemma 14.31 

A brief outline of the paper: We introduce some notation and state the 
main theorems in the next section. Section [3] is devoted to some auxiliary 
results including the L2-estimates. In Section |4] we prove some Holder esti- 
mates. We obtain the interior mean oscillation estimates in Section [5] and 
prove the solvability of systems in the whole space, i.e. Theorem [221 in Sec- 
tion m In Section [71 we establish several boundary estimates including the 
boundary mean oscillation estimates and prove the solvability of systems on 
a half space, i.e. Theorem 12.41 Finally we deal with parabolic and elliptic 
systems on a Reifenberg flat domain in Section [8l 

2. Main results 

Before we state our assumptions and main theorems, we introduce some 
necessary notation. By we mean a d-dimensional Euclidean space, a point 
in M*^ is denoted by x = (xi, . . . , x^) = {xi,x'), and {ej}j^^ is the standard 
basis of M''. Let M'^+i := M x M'^ = {{t, x) :t£R,x e R'^}. Throughout the 
paper, indicates an open set in R'^ and $7^ := (-oo,r) x C 1^"^+^ For 
vectors € C", we denote 

n 

(e,r?) = 

1=1 

For a function / defined on a subset P in M'^'^^, we set 

{f)v = 7^ I f{t,x)dxdt= / f{t,x)dxdt, 

\^\ JV JV 

where \D\ is the d + 1-dimensional Lebesgue measure of T>. 

Throughout the paper, we assume that the n x n complex-valued coeffi- 
cient matrices are measurable and bounded, and the leading coefficients 
\ce\ = 1/3 1 = m, satisfy an ellipticity condition. More precisely, we 
assume the following. 

(1) There exists a constant 6 € (0, 1) such that the leading coefficients 

|a| = = m, are bounded by and satisfy 

5iep< Yl '^Mt,x)c^,Co.)<5~'\^\^ (2.1) 

\a\ = \f}\—m 

for any {t,x) G M'^+^ and £, = {Ca)\a\=m, G C". Note that C can 

/m + d—l\ 

be considered as a vector in C <*-i where 

( d-1 j " ^ ^ " ml(d-iy. ■ 

I o I =m 

Here we use to denote the real part of /. 

(2) All the lower-order coefficients a^/j, \a\ ^ mov ^ m, are bounded 
by a constant K > 1. 



6 



H. DONG AND D. KIM 



Let 

Br{x) = {y € M"^ : |x - y| < r}, Bl{x') = {y' G M^^-^ : \x' - y'\ < r}, 

Q,{t, x) = {t- r^-, t) X Br{x), Q',{t, x') = {t- r^™, t) x 

On the leading coefficients we impose a very mild regularity assumption 
with a parameter 7 € (0,1), which will be specified later. To state this 
assumption, throughout the paper we write {aai3}\a\=\i3\=m ^ ^ whenever 
the n X n complex- valued matrices a^/? are measurable functions of yi G R 
only, {(iapl < 6~^, and {aai3}^a\=\i3\=m satisfies the ellipticity condition (|2.ip . 
For a linear map T from to M"^, we write T € O if T is of the form 

T{x) = px + ^, 

where p is a d x d orthogonal matrix and ^ € M'^. 

Assumption 2.1 (7). There is a constant Rq S (0, 1] such that, for each 
parabolic cylinder Q := (to — r^™,to) x Br{xo) with r < Rq, one can find 
Tq £ O and coefficient matrices {0.^(3} \a\=\/3\=m. ^ -'^ satisfying 

sup / \aai3{t,x) -aai3iyi)\dxdt <j\Q\, (2.2) 

\a\ = \l3\=mJQ 

where y = Tq{x). 

Let us introduce some function spaces utilized throughout the paper. We 
define the solution spaces ?{™((S', T) x Vl) as follows. For a domain Q. in M'^, 
we set 

]H;'"((S,T) X 0) = {/ : / = ^ G i^p((5,T) X 0)}, 

|o|<m 

II/IIh-'"((5,t)xC) = ™f { X] ll/"llip{(5,T)xn) : / = X] ^"/a}- 

|a|<m |o|<m 

Then 
?^^((S,T) xf^) 

= {n : txt e ]H-'"((5,T) x 0),Z)"n G Lp{{S,T) xn),0< \a\ < in}, 

\\u\\n^aS,T)xQ) = ll^tllH-™({S,T)xQ) + X] ll-^"^llip{(5,T)xn)- 

|a|<m 

In this paper u G C{^(X') means that u is infinitely differentiable on D, 
where 2? is a subset of either 1^"^+^ or R'^. As usual, C^{V) means the 
collection of infinitely differentiable functions with compact support (s "D. 
We define C^{[S,T] x 0) to be the collection of infinitely differentiable 
functions (^(t, x) defined on [S, T] x such that, for each t G [S, T], (f){t, ■) G 
C^{VL). The reader understands that if either S or/and T is infinity, then the 
closed interval should be replaced by a half-open or open interval. Finally, 
we set iq,\{S, T) x f)) to be the closure of C^{[S, T] x n) in ?^;"((S, T)xn). 
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Now we state the main result concerning parabolic systems in divergence 
form defined in the whole space. 

Theorem 2.2. Let n = R'^, p £ (l,oo), T G (-00, +00], and = 
{f^,...,f^Y^ G Lp{Qt), < m. Then there exists a constant 7 = 
^{d,n,m,p,5) such that, under Assumption \2. 1\ h). the following hold true, 
(i) For any u G 'H'^{VLt) satisfying 

ut + {-irCu+\u= ^"/a (2-3) 

\a\<m 

in VLt, "we have 

X^-^\\D^u\\L^(^nr)<N Y >^^\\fJLAnT)^ (2-4) 

a|<m |q|<??i 

provided that A > Aq, where N and Aq > depend only on d, n, m, p, 5, K 
and Rq . 

(a) For any A > Aq, there exists a unique u G T-L^{il.T) satisfying (j2.3p . 
(Hi) If all the lower-order coefficients of C are zero and the leading coeffi- 
cients are measurable functions of xi G M only, then one can take Aq = 0. 

Our next result is about the Dirichlet problem on a half space. For this, 
we impose the following assumption, where the parameter 7 G (0, 1 /4) is to 
be determined later. Set = {{xi,x') G M"' : xi > 0} and M^f+^ = M x R^. 

Assumption 2.3 (7). There is a constant Rq G (0, 1] such that the following 
holds with Q := {to - r'^"^,to) X Br{xo). 

i) For any x G M^J., t G M and any r G (0, min{i?o, dist(x, 90)}] so that 
Q C M^J.'^'^, one can find 7q G O and coefficient matrices {aai3}\a\=\i3\=m ^ ^ 
satisfying (j2.2p . 

ii) For any x G dW^, t G M and any r G {0,R(j], one can find Tq G 
O satisfying pu > cos(7/2) and coefficient matrices {aai3}\a\=\[3\=m ^ ^ 
satisfying (|2.2p . 

With a sufficiently small 7, the condition pn > cos(7/2) means that at 
any boundary point the yi-direction is sufficiently close to the xi-direction, 
i.e., the normal direction of the boundary. 

Theorem 2.4. Let Q. = R'^, p G (1, 00), T G (-00, +00] and 

/a = Lp(Ot), \a\<m. 

Then there exists a constant 7 = j{d, n, m,p, 6) such that, under Assumption 
\2.3\ (■y), the following hold true, 
(i) For any u G 'H'^{VLt) satisfying 

ut + {-lTCu+\u= Y ^"/a (2.5) 

|a|<m 
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in VLt, we have 

|Q:|<m |Q|<m 

provided that A > Aq, where N and Aq > depend only on d, n, m, p, 5, K 
and Rq. 

(a) For any A > Aq, there exists a unique u E T-U^IVLt) satisfying (j2.5p . 
(Hi) If all the lower-order coefficients of C are zero and the leading coeffi- 
cients are measurable functions of xi G M only, then one can take Aq = 0. 



Solutions of (j2.5p or (|8.2p below are understood in the weak sense: we 
say u € ^^{{S.T) x 0) satisfies ([23]) in {S,T) x n ii 

rT 

(j)t • udx dt 



[ cP{T, •) • uiT, ■)dx- f [ 
Jn Js Jq 

+ Y [ [ ((-l)™^'"'^"'^-aa/3^'^M+A0-u) (ixdt 
~i Js Jn Jn 



for any test function </> = 02, . . . , 0") G C^{[S, T] x 0). li S = -oo or 
T = oo, we assume (j){—oo, •) = or (j){oo, •) = 0, respectively. 

Remark 2.5. The ellipticity condition (j2.ip can be relaxed to a weaker 
condition 

d 

5^|eme,|'< Y 3ftK/3(*'^)e/3,ea) <<5~'l^l'- (2.6) 

j=l \a\ = \l3\=m 

For instance, when d = m = 2 the operator C = Df + Z)| satisfies (j2.6p 
with (5 = 1, but is not elliptic in the sense of (|2.ip . However, the condition 
(|2.ip has the advantage that it is invariant under orthogonal transformations 
of the coordinates. We claim that any operator £. satisfying (|2.6p can be 
rewritten into another divergence form operator which satisfies (j2.ip with a 
possibly different 6. In the above example, one way is to write 

Dfu + D^u = Dfu + D^u - DliD^u) + L»i2(Di2u). 

The symbol of the right-hand side is 

^{2,0) + ^{0,2) - ?{2,0)?(0,2) + 

which obviously satisfies (12. ip with 5 = 1/2. 

The claim is a simple consequence of the following observation. We only 
consider the case d = 2. The general case follows from an induction using 



HIGHER ORDER ELLIPTIC AND PARABOLIC SYSTEMS 9 

linear interpolations to cover the convex hull of d vertices. For simplicity, 
we also assume n = 1 and everything is real. First, it is easy to check that 

m m—1 

X] ^U,m-j) ~ X C{j-l,m-j+l)C{j+l,m-j-l) 
j=0 j=l 
m m—1 

- ^ ^{j,m-j) ^ ["^ ?(j-l,m-j+l) + ^ ) 

j=0 j=l 

m—1 m 

- X ?0-l,m-i+l)?{i+l,m-i-l) > ~^^ij,m-j)- 
j=l j=0 

Therefore, there exist e = e(m) > and 6i = 5i{m) > such that 

2 m-l 

Then suppose C satisfies ()2.6p . Using the fact that 

Jjih-rn-j) jj{j,m-j) _ jj(j-l,m-j+l)-Q{j+l,m-j-l)^ 

we then rewrite £ as 

m—1 

£ _|_ £ ^ (^f J)(hm-j) jj{j,m-j) _ 2j^^(j-l,m-j+l)^{j+l,m-i-l)^ ^ 
J=l 

which satisfies (|2.ip with (5(5i in place of 5. This completes the proof of 
the claim. Note that the leading coefficients of the new operator satisfy the 
same regularity assumption as those of C. Therefore, the results of our main 
theorems still hold true under the condition (12.60. 



3. Some auxiliary estimates 
In this section we consider operators without lower-order terms. Denote 

\a\ = \j3\=m 

3.1. L2-estimates. The first result is the classical L2-estimate for para- 
bolic operators in divergence form with measurable coefficients. We give a 
sketched proof for the sake of completeness. 

Theorem 3.1. Let T G (-oo, oo] and Q = W'-. There exists N = N{d, m, n, 6) 
such that, for any A > 0, 

X \'-^^\\D^u\\L.^^nr)<N E >^^\\fa\\L,inr), (3-1) 
|a|<m |a|<m 
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provided that u € W^IVLt), fa € L2(^^t)> |a| < "m, and 

Ut + {-irCou + Xu= D^fa (3.2) 

\a\<m 

in Qt- Furthermore, for any A > and f^ G L2{0,t), |a| < m, there exists 
a unique solution u € ^{^{^t) to the equation (j3.2p . 

Proof. By the method of continuity and a standard density argument, it 
suffices to prove the estimate ()3.ip for u € Cq°((— oo,T] x fi). From the 
equation, it follows that 

/ \{u,tLt) + {D°'u,aai3D^u) + \\u\'^ dxdt 

= (-1)'"' / iD''uJa)dxdt. (3.3) 

\a\<m -^^T 



By the uniform ellipticity (|2.ip . we get 

5/ \D"'u\'^dxdt < [ •^{aapD^u,D°'u)dxdt. 



We also have 



/ ^{u,ut)dxdt = - [ \u\'^{T,x)dx >0. 
Hence, for any e > 0, 
sf \D'^u\^ dxdt + \ I \u\'^ dx dt < J2 i-^)^''^ [ ^{D'^u,f^)dxdt 

^ ^ \(y.\<m ^ 

Em—\a\ / o -1 X — ^ m—\a.\ I ^ 

/ \D''u\^dxdt + Ne-^ ^ A — — / \fj^dxdt. 

I ^ ^ 11^ 

a|<m |a|<m 

To finish the proof, it suffices to use interpolation inequalities and choose e 
sufficiently small depending on (5, d, m and n. □ 

We also have the following local L2-estimate, where for a later use we 
consider only a simple case that the right-hand side of (|3.4p is zero. 



Lemma 3.2. Let^ <r < R<oo. Assume u € C'j^^iW-'^^) and 

ut + {-irCou=Q (3.4) 

in Qr. Then there exists a constant N = N{d,m,n,6) such that for j = 
12... m 

\\D^u\\L,^Q^^<N{R-r)-^\\uU,iQ^y (3.5) 

Proof. First we consider the case j = m. Set 

J^R-r , 
ro = r, rk = r + 2_^ , A: = l,2,..., 

1=1 



HIGHER ORDER ELLIPTIC AND PARABOLIC SYSTEMS 11 

.. = ^^^, ^ = 0,1,2,.... 



We choose nonnegative real-valued Ckit^x) G C^(M ) such that 
J 1 on Qr^, 

|0 on IR'^+^(-4-,4-)xS,,, 




and 



\iC^)^\<^J]^^^ l^'^^-l - ^(ii^T)!' ' = 0,l,...,m. (3.6) 
By applying ttC| as a test function to the system p.4p we get 

/ {uCl,ut)dxdt+ I {D°'{uCl),ao,i3D^u)dxdt = {). (3.7) 



Note that 



/ {u(l,Ut)dxdt= / |i«^fcp(0, x) — / {ut,uC,1) dx dt 
JQr Jbr Jqr 

2Ck{Ck)t\u\'^ dx dt, 



Qr 



which shows that 



^[ {iKlut)dxdt = l [ (0,x)dx- /" CkiCk)t\u\^dxdt. (3.8) 

■JQr ^ J Br JQn 

On the other hand, by Leibniz's rule 

/ (D^{'u/:l),a^^D''u) dxdt= [ (CkD'^u^a^^CkD^u) dxdt 
JQr ^ ^ JQr^ ' 

+ [ Yl Ca,,a2{D"'Ci){D'''u,aa^D^u)dxdt -.= 11+12, (3.9) 

JQr ai+a2=a 
\a2\<m 

where ai,a2 are multi-indices, and Cqj^qj corresponding appropriate 
constants. By the ellipticity condition (j2.ip . it follows that 

5?(/i) = / ^UkD'^u^a^pCkD^u) dxdt>5 f \CkD'^u\^ dx dt. (3.10) 

Here, we recall Mq = (—00,0) x W^. To estimate /2, we first see that 
and u in the integrals can be replaced by Mq and wCfc+i) respectively. Then 
using ()3.6p . we have 

™-l 2(m-l)k r 

\h\<NY, E / Ji?"^(Ma+i)||Z?^«fe+i)|dxdt:=/3. 



«=0 |a2|=«,|/3|=m 

Set 



b = ii«iiL(qh)- 
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Combining (13. 7p . ()3.8p . ()3.9p . and ()3.10p as well as using the inequality for 
(Cfc)i in p.6p . we obtain 

p r)2mk 
^ y^JC^-P'"^P'^^^^<^ (^_,)2>n B + ^3- (3.11) 

To estimate Is, using Young's inequality we observe that, for each < / < 
m — 1, 

\D\uCk+i)\\D"^iuCk+i)\dxdt 



where e > is an arbitrary real number. Furthermore, for / = l,...,m— 1, 
by interpolation inequalities 

/ \D'{uCk+i)\^dxdt<eo [ \D"'{uCk+i)\^dxdt + Nef^B, 

where we set Eq = 4e'^{R — 7')2{m-/)22(«-m)fc_ Combining the above two 
inequalities with (j3.1ip implies 

d[ \CkD"'u\'^dxdt<e [ \D"'{uCk+i)\'^ dx dt 

+ m \ 

— 771 I 

{R-r 



+ Ar 1+^ei^ B. (3.12) 



1=0 



Now we set 

^ \D'^{uCk)?dxdt. 



To estimate A^, we use (|3.6p and interpolation inequalities to get 



22mfc 



\uD^Ck?dxdt= / \uD^Ck?dxdt<N- .„ . 
and, for 1 < Z < m — 1, 

/ \D'^-'CkD'u\^dxdt= [ \D^-'CkD'{uCk+i)fdxdt 



^^ (^_,)2(^-0 ljD\^^.^)?dxdt 



<e \D"'{y^,^,)\^dxdt + Ne—— — B. 

By Leibniz's rule, we estimate Afc by 

Afc<eAfc+i + iV / iCfcD^txpdxdt + iV— ^ Ve^B. 

(^ii: - r) 



1=0 
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™-l , , , , , 92mA: 



This combined with (I3.12p shows 

m — 1 

We multiply both sides of the above inequality by and sum over k to 
obtain 

oo oo m— 1 oo 

k=0 k=l 1=0 k=0 

Choose e = 2~^™~^ and observe that X^^Lo^^^fc < cxd. Then the above 
inequality gives 

Ao < iV(i?-r)^2mg^ 

which clearly implies the desired inequality ()3.5p when j = m. 
Since 

Ao= / |D™«o)Pc^2;dt, 
the proof above shows that 

[ \D"'{uCo)\^dxdt<N{d,m,n,5){R-r)~^"' [ \u\^dxdt. 

Then for < j < m, by interpolation inequalities as well as the above 
inequality, 

/ \D^u\^dxdt< I \D^ {uCo)\^ dx dt < {R - r)-'^^ I \uCo\^dxdt 

+ N{R-rf^"'-^^ [ \D"'{uCo)\'^dxdt<N{R-ry^^ [ \u\^dxdt, 
Jr^ Jqh 

where N = N{d, m, n, 6). The lemma is proved. □ 

We are going to use the following Poincare type inequality, which gener- 
alizes a result in [271 Sect. 3] 

Lemma 3.3. Let p G [l,oo), R e (0, oo), u G C^^{W^+'^). Suppose that u 
satisfies 

uj + (-ir£ix = (3.13) 

in Qr. Let P = P{x) he the vector-valued polynomial of order m — 1 such 
that 

{D''P)q^ = {D''u)q^, k = 0,l,...,m-l, 
and let V = u — P. Then for each k = 0,1, . . . ,m — 1, we have 

WD'vU^^Qn) < iVi?"-'^P"«||L,{Q,), (3.14) 
where N = N{d, m, n, 6,p) > 0. 
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Proof. By a simple scaling, without loss of generality we may assume that 
R = 1. Take a function € Cq°{Bi) with unit integral. For any k = 
0, 1, . . . , m - 1 and t G (-1, 0), let 

9kii) = [ Ciy)D''vit,y)dy. 



Then for any t € (—1,0), by Holder's inequality and Poincare's inequality, 
we have 

/ \D''v{t,x)-gf,{t)\Pdx= f \f {D''v{t,x)-D''v{t,y))C{y)dy\''dx 
<N [ [ \D^v{t,x)-D^v{t,y)\Pdydx<N f \D^+^v{t,x)\P dx. 

JBi Jbi JBi 

(3.15) 

Now let Cfc = Qkit) dt be a constant vector. Since 

D'^vdxdt = 0, 

Qi 

by the triangle inequality, ()3.15|) . and Poincare's inequality, we get 

< N\\D''v-Ck\\L,iQ,) ^ ^ll^'^-9fcllLp(Qi)+^ll9fc-Cfe||L,(Qi) 

< N\\D^+'vU^(^Q^)+N\\dtgk\\LM-i,o))- (3-16) 
By the definition of Qj^, (13.130 and integration by parts, 

dt9kit)= I C{y)D''dtv{t,y)dy= [ ay)D''dtu{t,y) dy 

JBi JBi 



\m+l 



Bi 



ay)D^Cuit,y)dy 



= (_l)fc+i / (D''D''C)iy)aa^Dpuit,y)dy. 
JBi 

Thus, by Holder's inequality, for any t G (—1,0), 

\dtgdt)\<N\\D"^u{t,.)\\LAB,y (3.17) 

Notice that D'"^v = D'^u. Combining (f3T6]) and (fSTTTl) yields (|3TT4l) by an 
induction on k. □ 

Now we prove an estimate for D'^u when a = (ai, . . . , a^) satisfies \a\ > 
m and ai < m. 

Corollary 3.4. Let < r < R < oo and aa/s = CLafsixi), \a\ = = m. 
Assume that u G Cf^^iW^'^^) satisfies (|3.4p in Qji. Then for any multi-index 
a = (ai, . . . , ad) such that \a\ > m and a\ < m, we have 

where N = N{d, m, n, 6, R, r, a) . 
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Proof. Since D^',u also satisfies (|3.4p . by applying Lemma [3.21 repeatedly, 
we obtain 

||I)"tt||L.(Q.) < iV E \\D''^\\L2iQn) (3-18) 

k<m 

for any a with ai < m. Now let P = P(x) be the vector-valued polynomial 
of order m — 1 such that 

D'^p) =(d^u) , = 0,1,... ,m- 1, 

and set v = u — P. Then we have 

D^v] =0, A: = 0,1,... ,m- 1. 



Since v satisfies (13. 4p in Qr, by (I3.18P with in place of u we obtain 



k<m 



for any a satisfying |q| > m and ai < m, where the second inequality is 
due to Lemma 13.31 □ 



Lemma 3.5. Let < r < R < oo and a^p = aap{xi), \a\ = |/3| = m. 
Assume that u € C'^^{M.'^'^^) satisfies ()3.4p in Qr. Then for any integers 
i > 1, J > 0, and any multi-index a = (ai, . . . ,arf) such that \a\ > m and 

ai < m, we have 

+ mo^uu^i^Q^^ < iviiD'"t.iu,(Q^), (3.19) 

where N = N{d, m, n, 5, R, r, a, i, j) . 

Proof. Since iLt satisfies ()3.4p in Q^, by Lemma 13.21 

II^™'*^IIl2(q,,) < ^^IKIIl2(q,2)' r<ri<r2<R. 
Using this inequality, Corollarv 13. 4| and the fact that dlD^',u also satisfies 



(j3.4p in Qfi, we see that, in order to prove ()3.19p . it is enough to show 

WMi.iQr) < N{d,m,n,6,R,r)\\D^u\\L^^Q^), (3.20) 
where R may be a smaller one than that in ()3.19p . Take r^, s^, and Cfc from 



the proof of Lemma 13.21 Set 



Q^''^ = Qr„ Q^''^ = Qs„ k = 0,1,2,.... 



Also set 

A — lU. I|2 R _ II n'"--i.l|2 

We now apply tkCl *° (|3.4p as a test function to get 



Ic 



UtQk?dxdt = -j {D°'{utQl),aapD^u)dxdt. (3.21) 
Qr JQr 
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To estimate the terms in the right hand side, we first observe that, thanks 
to the fact that ut also satisfies (13. 4p in Qr, by applying p.Sp with Sk and 
rfc_|_i in place of r and R, respectively, 

/ \D^Ut\^ dxdt<N{rk+i- SkY"^^ I \ut\^dxdt 

= ^Tn ToT / Itkl"^ dxdt, I = 0,1, . . . ,m. 

{R - ry^ jQ(k+i} 

Hence, for each / = 0, 1, . . . , m, by Young's inequality 

\D''^-^{Cl)\\D\\\D'''"u\ dxdt 



2{m-l)k / f If 

<- — jrieo \D^ ut\^ dx dt + — \D'^u\^dxdt 

<£ IzLtl^ dxdt + N-- --— / iD'^ul^dxdt, 

where e > is an arbitrary real number. This along with (|3.2ip and Leibniz's 
rule implies that 

Finally, we prove (j3.20p following the same argument as in the proof of 
Lemma 13.21 □ 



3.2. Maximal and sharp functions. We recall the maximal function the- 
orem and the Fefferman-Stein theorem. Let 

Q = {Q.^{X) : X = {t,x) G]R'^+\r G (0,oo)}. 

For a function g defined in M*^"*"^, the (parabolic) maximal and sharp function 
of g are given by 

Mg{t,x)= sup / \g{s,y)\dyds, 
QeQ,{t,x)eQ JQ 

g*{t,x)= sup / \g{s,y) - {g)Q\dyds. 



It is well known that 

ll5'llLp(M'i+i) < ^h^WLpiRd+i), ||-^5||Lp(IR'i+i) < ^I|5'IIlp(R'«+1)> 

if 5 € Lp(M'^+^), where 1 < p < oo and = N{d,p). Indeed, the first of the 
inequalities above is due to the Fefferman-Stein theorem on sharp functions 
and the second one to the Hardy-Littlewood maximal function theorem (this 
inequality also holds trivially when p = oo). 

Theorem 13.61 below is from |29j and can be considered as a generalized 
version of the Fefferman-Stein Theorem. To state the theorem, let 

Q = {C/(io,n, ■ ■ .,id),io,k,- . . ,id G Z}, / S Z 
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be the collection of partitions given by parabolic dyadic cubes in M.'^^^ 

[io2-2-', (io + 1)2-2-') ^ [.^2-', ih + 1)2-') X ... X [z,2-', (i^ + 1)2-')- 

Theorem 3.6. Let p £ (l,oo), and U,V,F e Lijod'^^'^^^)- A ssume that we 
have \U\ < V and, for each I G Z and C € C;, there exists a measurable 
function on C such that \U\ < U'-' <V on C and 

[ \U^ - (U^) \dxdt < [ F{t,x)dxdt. 
Jc Jc 

1^ < iV(rf,p)||F||i^(Kd+i)||F||^^J(j,+i), 



Then 



provided that F,V Lp 



4. Interior Holder estimates 

By using the L2 estimates obtained in Section [3.11 in this section we shall 
derive interior Holder estimates of derivatives of u. As usual, for /_f € (0, 1) 
and a function u defined on P C M*^"^^, we denote 



\u{t,x) - u{s,y)\ 



(t,x),{s,y)£V |t ~ Sr' -\- \X — y\f^ 



Lemma 4.1. Let aa/s = aa/sixi). Assume that u G C^^M.'^'^^) satisfies 
(j3.4p in Q2. Then for any a = (ai, . . . , a^) satisfying \a\ = m and ai < m 

we have 

where N = N{d, m, n, 6) > 0. 

Proof. Thanks to the well-known interpolation inequality, it is sufficient to 
estimate [D°'u](~<i/2(^Q^y The proof is based on a convenient form of the 
Sobolev inequality. By the triangle inequality, we have 

\D'^u{t,x)-D'^u{s,y)\ ^ ^ , ^ 
sup — —J- j , < 1 + J, 

{t,x),{s,y)eQi |t - s|V4 + |a; _ 

{t,x)^{s,y) 



where 



^ _ \D'^u{t,xi,x') - D'^u{t,yi,x') 



I 11/2 ' 

xi,yi£{-l,l),xiytyi Fl-yi|' 
{t,x')€Q[ 

\D"u{t,y,,x')-D"u{s,yi,y') 



J := sup 



yimi) |t-s|V4 + |a,'-y'|l/2 

it,x'),{s,y')GQ[,it,x')^{s,y') 
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Estimate of I: By the Sobolev embedding theorem D°'u(t, 
function of xi G (—1, 1), satisfies 

\D''u{t,xi,x') - D°'u{t,yi,x')\ , 

^,"-1,1) N^^^F 

(4.1) 

On the other hand, there exists a positive integer k such that D°'u{t,xi,x') 
and DiD°'u{t,xi,x'), as functions of {t,x') G Q'^, satisfy 

sup (jZ)°M(i,xi,x')| + \DiD'^u{t,xi,x')\) 

<iV||Z)"t.(.,Xi,.)|liy.(Q;)+iV||I?lI)"«(-,Xi,-)|| H/|(Q'i)- 
The inequahty above imphes that, for ah {t,x') G Q'l, 

j \D'^u{t,xi,x')\^ dxi + j \DiD'^u{t,xi,x')\^ dxi 

i+\P\<m+l+k 
I3i<m 

This combined with (14.11) shows that 



I<N Yl \\9lD"uh^(^Q^) < iV||Z)-«|U,(Q,), 

j+|/3|<m+l+fc 

where the last inequahty is due to Lemma |3.5[ 

Estimate of J: By the Sobolev embedding theorem, we find a positive 
integer k such that D°'u{t,yi,x'), as a function of {t,x') G Q'l, satisfies 

(4.2) 

For each i,j such that i + j < k, dlD-^^,D°'u(t,yi,x'), as a function of yi G 
(—1,1), satisfies 

sup \diDi,D''u{t,yi,x')\ 
?/ie(-i,i) 

This together with (14. 2p and Lemma 13.51 gives 



i+\l3\<m+l+k 
/3i<m 

This completes the proof of the lemma. □ 
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For A > 0, let 



U= ^ A5-^|D"u|, U' = J2 A^-5^|D"ii|. (4.3) 

|o|<m \a\<'m,ai<m 

Corollary 4.2. Let 0^13 = aa/3(xi) and A > 0. Assume that u G C'^^{M.'^~^^) 
satisfies 



ut + (-l)™Xoit+ Aii = 



(4.4) 



in Q2. Then we have 

where N = N(d, m, n, 6) > 0. 

Proof. The case when A = follows from Lemma 14. li To deal with the 
case A > 0, we follow an idea by S. Agmon. Let r]{y) = cos(A^/^^"*^y) + 
sin(A^/^^"^^y) so that rj satisfies 

= (-l)-Ar?, r?(0) = l, |Z)^r?(0)|=A^/(2-) j = l,2,.... 



Let z = {x, y) be a point in M^"^^, where x E W^, y E IK, and u{t, z) and Qr 
be given by 

u(t, z) = u{t, X, y) = u{t, x)r]{y), Qr = {-r'^'^, 0) x {\z\ < r : z e M''+^}. 
Since ii satisfies, in Q2, 

iit, + {-ircou+{-irDi"^{u) = o, 

by Lemma 14.11 applied to ii we have 

^'HlcV.(QO - mni,n,6)\\DTu\\,^^Q^^ (4.5) 

for any /3 = (/3i, . . . , P^+i) satisfying \f3\ = m and /^i < m. Notice that for 
any a = (ai, . . . , a^) satisfying \a\ < m and ai < m, 



1 \a\ 
A2 2,71 ||D"lt||„i 



Cl/2{Qi) 



< N 



Cl/2(Qi) 



^ = (ai, . . . ,ad,m - \a\ 



and D^u is a linear combination of 



A2 2m cos(A2my)L)^u, A2 2m sin( A 27n y ) u, /c = 0, 1, . . . , m. 



Thus the right-hand side of (j4.5p is less than the right-hand side of the 
inequality in the lemma. The lemma is proved. □ 

Let a = mei = {m,0, ... ,0). In the remaining part of this section, we 
shall establish a Holder estimate of 



m=rn 

We make use of the following elementary lemma. 



(4.6) 
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Lemma 4.3. Let r G (0,oo), k > 1 be an integer, p G [l,oo], and u £ 
Lp{[0,r]). Assume that D^u = fo + Dfi + ... + D^^^ft _i in (0, r), where 
fj G Li{[0,r]),j = 0,...,k-2 and fk^ie Lp([0,r]). Then Du G Lp([0,r]) 
and 

k-2 
j=0 

where N = N{k, r) > 0. 

Proof. Thanks to scaling, we may assume r = 1. For any function / G 
Li([0, 1]), we define its anti-derivative Zf : [0, 1] ^ M as X/(x) = f{t) dt. 
It is easily seen that D'^u = D^^^fk^i, where 

fk-l = fk-l + + ^"^fk-S + • • • + I'^'^fo, 

and 

k-2 

ll/fc-l||Lp([0,l]) < ^ll/fc-lllLp([0,l]) + ll/7'llii([0,l])- 

j=0 

Therefore, without loss of generality we may assume fj = for j = 0, . . . ,k— 
2. Under this assumption, we have for some constant Cj,j = 0,l,...,k—1, 

u{x) = {lfk-i){x) + Co + 2cix + . . . + kck-ix^-^. (4.8) 
We claim that 

|c,-| < A^||n||i^([o,i])+iV||/fc_i||i,([o,i]), i = 0,l,...,/c-l, (4.9) 

which immediately yields (j4.7p . To prove the claim, we integrate both sides 
of (gSI) on [0,i/A;], j = 1, 2, . . . , A; to get 

H/k f-j/k 

/ u{x)dx= / {Ifk-i){x)dx + coij/k)+ci{j/k)^ + ... + Ck-i{j/kf. 
Jo Jo 

The claim (j4.9p then follows since the matrix [{j /ky]f j^-^ is nondegenerate 
and 

/ {u{x) - {Ifk-i){x))dx 
Jo 

< N\W\\li{[o,i]) +^ll/fc-i||Li([o,i]), J = 1,2, . . . , A;. 

□ 

Corollary 4.4. Let k > 1 be an integer, r G (0, oo), p G [l,c>o], P = 
[0,r]'^+i, and u{t,x) G Lp(D). Assume that D\u = /o + Di/i + ... + 
Di~^fk-i in Z), where fj G Lp{'D),j = 0, . . . ,k — 1. Then Diu G Lp(T>) and 

k-l 

\\DiuU^^r>)<N\\uU^tj,)+NY,\\fjKiv): 

j=0 

where N = N{d, k, r) > 0. 
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Proof. The corollary follows from Lemma 14.31 by first fixing {t,x') and then 
integrating with respect to {t,x'). □ 

Lemma 4.5. Let < r < R < oo and = aap{xi). Assume u € 
(j<xj^d+i\^ gdfigflgg (j3,4p in Qj^. Then, for any nonnegative integers i,j, 

WdlDieU^iQ,.) + \\dlDiD^QU,iQ^^ < N\\D"^uU^^Q^^, (4.10) 

where N = N{d, m, n, r, R, 5, > 0. 

Proof. Obviously, we have 

\\eh,iQ^) < NWD^uU^^Q^y (4.11) 

Thus as noted at the beginning of the proof of Lemma \'6.5\ it suffices to 
prove 

\\D,Qh^^Q^)<N\\D"^uh,^Q,^), (4.12) 
where R' = (r + ^)/2. From (|3.4p . in we have 

D'['e = i-ir+^ut- Yl D^iaa^Dpu) 

\a\=m=m 
ai<m 

\a\ = \(5\=m 

where a = (ai, a2, . . . , a^) = {ai,a'). Then the estimate (|4.12p follows 
from Corollary 14.41 with a covering argument and Lemma 13.51 The lemma 
is proved. □ 

The following Holder estimate is deduced from Lemma 14.51 in the same 
way as Lemma 14.11 and Corollary 14.21 are deduced from Lemma 13.51 

Lemma 4.6. Let a^^ = Oa^ixi) and A > 0. Assume that u € Cf^^ 
satisfies ()4.4p in Q2 . Then we have 

||S||C1/2(Q,) < iV||f/||L2(Q2)' 

where N = N{d, m, n, 6) > 0. 

Since the matrix [a^ia]i'j=i is positive definite, we obtain the following 
estimate by using Corollary 14.21 and Lemma 14.61 

Corollary 4.7. Let Oap = cLapixi) and A > 0. Assume that u G C;^(M'^"'"^) 
satisfies (14. 4p in Q2. Then we have 

||f^||L^(Qi) < A^||C^||l2(Q2)' 

where N = N{d, m, n, 6) > 0. 
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5. Estimates of mean oscillations 

Recall the definitions of U, U', and in (j4.3p and (j4.6|) . respectively. With 
the preparations in the previous section, we obtain the following estimates 
of mean oscillations of U' and 0. 

Lemma 5.1. Let r € (0, oo), k G [2,oo), A > 0, and aap = aapixi). 
Assume u € C'^^{W^'^^) satisfies 

Ut + {-l)'^CQU+\u = Q 

in Qkt- Then we have 

{W - {u')qA)q^ + (le - (e)Qj)Q^ < iVK"V2(f72)V2^, (5.1) 

where N = N{d, m, n, 5) > 0. 

Proof. By a scaling argument, we may assume r = 2/k. Then by Lemma 
T6l we have 

(I© - (0)qJ)q. < Nr^'mc^r^iQ,) < Nn~^l\u')]l\ 



The first term on the left-hand side of (|5.ip is estimated similarly by using 
Corollary □ 



For/, = (/i,...,/^r, we denote 



\a\<m 

Lemma 5.2. Let r G (0, oo), k G [4, oo), A > 0, G L2^ioc{'^'^^^), \a\ < m, 
and Oafs = aap{xi). Assume u G C'^^{R'^^'^) satisfies 

Ut + {-irCou+Xu= D'^L 

\a\<m 

in Qur- Then we have 

{\U' - (C/')qJ)q. + (le - (e)gJ)Q^. < NK''/\U')'^l+Nn"^+HF')l/l 

(5.2) 

N-^U <U' + e < NU, (5.3) 
where N = N{d, m, n, 5) > 0. 

Proof. The inequality (j5.3p follows from the definitions of J7, U' and as well 
as the fact that [aia]^j=i is positive definite. To prove (|5.2p . we adapt the 
idea in the proof of Theorem 7.1 in [28] and take into account the presence 
of A. We can certainly assume that u and have compact supports. In 
addition, we assume that Oq/j and are infinitely differentiable. If not, we 
take the standard mollifications and prove the estimate for the mollifications. 
Then we can pass to the limit because the constant in the estimate (15. 2p 
is independent of the regularity of a^/s and . 
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Take a C G such that 

C = 1 on C = outside (-(Kr)^'", {nrf"^) x B^r- 

By Theorem 13.11 for A > 0, there exists a unique solution nf {W^+^) to 
the equation 

Wt + {-irCow + Xw= D'^iCfa)- 

\a\<m 

Since all functions and coefficients involved are infinitely differentiable, by 
the classical parabolic theory, w is infinitely differentiable. The function 
V .= u — w is also infinitely differentiable and satisfies 

Vt + {-l)"'Cov + Xv = in 

We define V, V' , W, and W' in the same way as U and U' . Thus by Lemma 
[O(note that n/2>2) 

- (OqJ)q. + (le - (e)Qj)^^ < NK~yW)]^l^^, (5-4) 

where is defined in the same way as Q with u replaced by v, i.e. 

O := a^iiD^v, a = mei. 

\l3\=m 

Next we estimate w. By Theorem 13. II we have 

a|<r?i |a|<m 

In particular, 

{wX^<N.-^i[F-)Z, MQl<A^(i^^C- (5-5) 

Now we are ready to prove (j5.2p . From (j5.4p and (jS.Sp . and the fact that 
u= w + V, we bound the left-hand side of (j5.2p by 

d^' - iy')Qr\)Q. + (I© - + N{w')t 

which is less than the right-hand side of (j5.2p . □ 

Recall that the ellipticity condition (12. ip is invariant under any orthogonal 
transformation of the coordinates. 

Corollary 5.3. Let r G (0,00), k G [4, 00), A > 0, G L2,ioc(IR'^+^), 
I a I < ni, and a^p = 00/3(2/1), where p is a d x d orthogonal matrix and 
y = px. Assume u G C'^^iW^^'^) satisfies 

Ut + {-irCou + \u= D'^L (5.6) 

la|<r?i 
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in Q := Q^r- Then there exist a function depending on Q, and a constant 
N = N{d, m, n,5) > such that 

N~^U <U^ < NU, (5.7) 
Proof. Since u satisfies (|5.6p . we see that v{t,y) := u{t,p~^y) satisfies 

\a\<m 

in Q, where a^/j are the corresponding coefficients in the y-coordinates and 
y), \a\ = k, is a hnear combination of fa(t: P~^y), |o| = k. 
Set 

= V' + G, 

where V' and are defined as U' and in (|4.3p and (|4.6p . but with v and 
aa/3 in place of u and aajj, respectively. Then since the new operator also 
satisfies (|2.ip . the corollary follows from Lemma l5.2i □ 

6. Proof of Theorem 12.21 

In this section we finish the proof of Theorem 12.21 First we observe that 
by taking a sufficiently large Aq and using interpolation inequalities, we can 
move all the lower-order terms of Cu to the right-hand side. Thus, in the 
sequel, we assume all the lower-order coefficients of C are zero. Recall the 
definitions of O and A above Assumption 12.11 the constant Rq from 
Assumption 12.11 

Theorem 6.1. Let j G (0,1), A > 0, and T,a G (l,oo), l/r + l/a = 1. 
Assume u G C^^{W^^^) vanishes outside Q-yR^ and 

ut + {-irCu+Xu= 

\a\<m 

in QKr{Xo), where G L2joc{^'^~^^) ■ Then under Assumption \2.1\ (j), for 
each r G (0, oo), k > 4, and Xq := {tQ^xo) G M'^"'""'^, there exists a function 
depending on Q := QKr{Xo), such that we have (15. 7p and 

- (f/«)o,,A,,l)e,,;,,, < Nn-''^(uXlt,,, 

where N = N(d, 5, m, n, r) . 

Proof. Fix K > 4, r G (0,oo), and Xq G M'^"'"-^. First we consider the case 
when Kr < Rq. For Q = QuriXo), from Assumption (j2.ip (7), we find 
7q G O and {aap}\a\=\i3\=m ^ ^ satisfying (|2.2p . Then we see that u satisfies 

^ + (-l)"^D''{a^pD^u) + Xu= Y D'^L 

\a\ <m 
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where UajS = aa/sivi), y = 7q(x), and 

fa=fa + l|a|=m ^ ("l)" (Oa/? " ^a/?) 
|/3|=m 

Using Corollary 15.31 with a shift of the coordinates, there exists a function 
U'^ satisfying (j5.7p such that 

(If" - (C«)q,,x.,I)«,,,., < ""-"^(f^jil,..) + 

where = A^((i, m, n, 6) and F is defined by using in the same way as 
F. Observe that for |a| = m 

f \fa?dxdt<N [ \f\^dxdt + NI, (6.3) 

where 

1= \2 / \{aai3 — aai3)D^ u\'^ dx dt. 



|/3|=m 



By Holder's inequality, we have 

I < NjI'^'jI/^ (6.4) 

where 

Ji = ^ / \aap - aaj3\'^'' dxdt, J2 = \D™-u\'^'^ dx dt. 

Since nr < Rq, by Assumption 12.11 

Ji<N 1 |a„/3 - a„^| dx dt < N{Kr)'^+^"'j, 

where N depends only on d, m, n and 6. From the above estimates for Ji as 
well as the inequalities (|6.2p . (j6.3p . and (j6.4p . we conclude (|6.ip . 

In case kt > Rq, we take t/*^ = U. By the triangle inequality and Holder's 
inequality, the left-hand side of (|6.ip is less than 



< A'«'"*«(l««.)i'lTi,)(^")Q'iU) S A''<"'*«7'"-'(f/-)^'lU,, 
where = A^((i, m, n, 5). The theorem is proved. □ 

Corollary 6.2. Let 7 G (0,1), A > 0, and T,a e (l,oo), 1/r + 1/a = 1. 
Assume that u € Co"(M''+^) vanishes outside Q-^r^ and satisfies 

Ut + {-lTCu+\u= ^"/a 

\a\<m 
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in M'^"'"^, where € L2^ioc{^'^~^^) ■ Then under Assumption \2.1\ ('-/), for each 
/ S C S C;, and k > 4, there exists a function U'-'' depending on C such 
that (|5.7|) is satisfied and 

{\U^-{U^)c\)c<N{F,)c, (6.5) 

where N = N{d, 5, m, n, r) and 

F,{t,x) = K-'/\M{U')Y^' 



Proof. For each k > 4 and C Ci, let Qr(-'^o) be the smallest cylinder 
containing C. From Theorem 16.11 we fine with Q = Q^rl-^o)- Take 
= U^. Then by Theorem O as well as the facts that C C Qr{Xo) and 
the volumes of C and QriXo) are comparable, we have 

{\U<'-{U^)c\)c<N{d)I, 
where I is the right hand side of the inequality ()6.ip . Note that, for example, 

t2\ / A^/rr2^ 



for any X = {t,x) G C Thus / is less than a constant times Fi^{X) for any 
X € C, especially, it is less than a constant times {F^)fj. Hence we arrive 
at the inequality (j6.5p . This finishes the proof of the corollary. □ 



Theorem 6.3. Let p € (2,oo), A > 0, and G Lp(R'^+i). There exist 
positive constants 7 G (0,1) and N, depending only on d, 5,m,n,p, such 
that, for u G C^(M'^'^^) vanishing outside Q^Rq and satisfying 

ut + {-irCu + Xu= D'^fa, 

\a\<m 

we have 

l|t^llLp(Rd+i) < ^I|-^IIlp{R''+1)' 

where N = N{d,5,m,n,p). 

Proof. Let 7 > and k > 4 be constants to be specified below. Take a 
constant r such that p > 2t > 2. For each / G Z and C G C;, let U'^ be the 
function from Corollarv 16.21 We know that there exist positive Ni, i = 1,2, 
depending only on d, m, n, and 5, such that 

U < NiU^ < N2U =: V 

This along with Corollarv 16.21 and Theorem 13.61 implies that 

\\U\\l^ < N\\F,UJV\\1-' < N\\F,UjU\\l-\ 
Here we denote Lp = Lp(R.'^~^^). The above inequalities readily give 

\\U\\L<N\\F4r. (6.7) 
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From the definition of i^^ ()6.6I) and the Hardy-Littlewood maximal function 
theorem (recall that p > 2t > 2) it follows that 

Using the above inequality and (j6.7p . we have 

\\U\\l, < N^i~y^\\U\\L, + iVK-+i||F||L^ + 7VK'"+iy/(2-)||f/|l^^. 

It only remains to choose a sufficiently big k, then a sufficiently small 7 so 
that 

iV^-i/2 + ^r^-+iy/(2a) < 1/2. 

□ 

Proof of Theorem \2.2l Due to the duality argument and Theorem 13.11 it is 
enough to consider the case p > 2. We first prove the first two assertions 
for T = +00. In this case the estimate (j2.4p is proved using Theorem 16.31 
and the standard partition of unity argument. Then Assertion (ii) follows 
from the method of continuity and the existence of solutions to systems with 
simple coefficients, for instance, aa/s = ^aplrnxm- For T < cxo, we extend / 
to be zero for t >T, and then find a unique solution u € 'Hp'' {^00) of (j2.3p 
in Ooo) the existence of which is guaranteed by the argument above. This in 
turn also yields the existence of a solution of (|2.3p in satisfying (j2.4p . For 
the uniqueness, let u G 'H^{Qt) be a solution of (j2.3p with zero right-hand 
side in Qt- Take ii to be the even extension of u with respect to t = T. 
Then ii £ Ti^i^oo) satisfies (j2.3p in Q^o with the right-hand side vanishing 
for t < T. It is easily seen from the method of continuity that u = for 
t < T. Finally, Assertion (iii) is due to a standard scaling argument. □ 

7. Boundary estimates 

This section is devoted to the Dirichlet problem for (jl.lh . We shall follow 
the lines of Sections [3H5] to carry out the corresponding boundary estimates. 
As before, we denote Co to be the operator with zero lower-order coefficients. 
We introduce a few more notation. For any t € M, x G M'^ and r > 0, denote 

dist(a;, 9^2) = inf |x — y|, i^ri^) = ^ Br{x), 

Cr{t,x) = {t-r'^"^,t) X nr{x). 
B+{x) = Br{x)r\M.%, Qt{t,x) = {t-r^'^,t) x 

7.1. Boundary L2- estimates. Similar to Theorem 13.11 we have the fol- 
lowing L2-estimate on a half space or on a domain. 

Theorem 7.1. Let T £ (—00, 00] and Q be a half space or a domain, 
(i) There exists N = N{d,m,n,5) such that, for any A > 0, 

^ Al-^||Z)"«||i,(f,,) <N X^Wfah^iUr), (7.1) 

|Q|<m |a|<m 
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provided that u S 'H^{Q,t), fa ^ L2{0,t), \a\ < m, and in Q,t 

ut + {-irCon+Xu= D'^fa- (7.2) 

\a\<m 

(a) For any A > and G L2{^t), W\ < m, there exists a unique solution 
u € 'H2^{0,t) to the equation (j7.2p . In the case that Vl is a hounded domain, 
one can take A = 0. In this case, N depends on (7 as well. 

Proof. The proof is similar to that of Theorem [3TTJ The last assertion follows 
from the Poincare inequality. □ 

In the proof of Proposition I7.10[ we will extend Theorem 17.11 to allow 
a more general right-hand side. In the remaining part of this section, we 
consider systems only on a half space. 

First we prove the following boundary L2-estimate, the proof of which is 
almost the same as that of Lemma 13.21 and thus omitted. 



Lemma 7.2. Let < r < R < oo. Assume u € C^^iR^^^) satisfies 

D^u = ondR'l, k = 0,1, . . . ,m - 1, (7.3) 

and 

ut + i-irCou = (7.4) 
in Q^. Then there exists a constant N = N{d,m,n,6) such that for j = 



1 m, 



\D'u\\^^^Q^^<N{R-r) -|i"iii,(Qj). 



Corollary 7.3. Let < r < R < oo and = CLctfiixi), \a\ = = m. 

Assume that u £ Cf^^{R'^'^) satisfies ([73]) and ([731) in Q^. Then for any 
integers i > 1, i > 0, and any a satisfying ai < m, we have 

where N = N{d, m, n, 5, R, r, a, i,j). 

Proof. Since D^,u also satisfies (|7.3p and (|7.4p . by applying Lemma 17.21 
repeatedly, we obtain for any a satisfying ai < m, 

where R' = {r + R)/2. From this inequality and the boundary version of the 
Poincare inequality along with the zero boundary condition (j7.3p it follows 
that 

The get the desired estimate, it suffices to use the argument in Lemma 
[331 □ 

Recall that 

G := a^pD^u, a = mei. 

|/3|=m 
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Lemma 7.4. Let < r < R < oo and = aQ^(xi). Assume u G 

C'^^iM.'^j^^) satisfies (|7.3p and ()7.4p in Q^. Then, for any nonnegative inte- 
gers i and j, 

where N = N{d, m, n, r, R, 6, i,j) > 0. 
Proof. Obviously, we have 

I|0|Il.(q+)<^^P"«IIl.(q+)- (7.5) 
As before, it suffices to prove that, for R' = (r + R)/2, 

< )■ (7.6) 

From (j7.4p . in we have 

\a\ = \l3\=m 
ai<m 

= (_l)-+i^_ ^ Dl^{a^pD^:D^u). 

\a\ = \(3\=m 
oi<m 

Then the estimate (j7.6p follows from Corollarv 14.41 with a covering argument 
and Corollarv 17.31 The lemma is proved. □ 

7.2. Boundary Holder estimates and Hardy's inequality. In the same 
fashion as in Section [H we obtain 



Lemma 7.5. Let = aap{xi). Assume that u G C~(M:|+^) satisfies 
()7.4p in Q2 and (j7.3p . Then for any a satisfying \a\ = m and ai < m, we 
have 

where N = N{d, m, n, 6) > 0. 

Proof. We follow the proof of Lemma 14.11 by using Corollary 17.31 instead of 
Lemma 13.51 □ 



Corollary 7.6. Let a^p = aap{xi) and A > 0. Assume that u G C^^iE-i^^) 
satisfies (j7.3p and 

ut + (-l)'"£o«+ A« = 

in Q2. Then we have 

\\U'\\ci/2^Q+) < N\\U\\^^^Q+y 

where N = N{d,m,n,6) > 0, and U,U' are defined as in (|4.3p . 

Proof. The corollary follows from Lemma 17.51 and Agmon's idea as in the 
proof of Corollary 14.21 □ 

As an analogy of Lemma 14.61 we obtain 
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Lemma 7.7. Under the conditions of Corollary \7.6[ we have 

II©IIC1/2(Q+) < ^\\U\\l2^q+), 

where N = N{d, m, n, 6) > 0. 

Similar to Corollary 14.71 Lemma 17.71 and Corollary 17.61 imply 
Corollary 7.8. Under the conditions of Corollary \ 7. 6[ we have 

where N = N(d, m, n, 6) > 0. 

We will use the following Hardy type inequality. 

Lemma 7.9. Let R € (0, oo), p G (l,oo], m be a positive integer, and 
f G C^{[0,R]). Suppose that /(O) = Df{0) = ... = D'^-^fiO). Then we 
have 

||x-'=Z)™"V(^)||^^([o,H]) <iV||Z)™/||^^([o,jjj), k = l,...,m, (7.7) 

where N = N{p,k) > is a constant. Furthermore, for any function g G 
Lg{[0, R]), q = p/{p - I), and r, G C£^((0, R]) satisfying 

\D''ri{x)\ <A'x-^ x G (0,i?], A; = 0, K > 0, 

we have 

ll^'"(^/)llL,([0,iJ])<A^II^'"/llM 
\\9D'^{vf)\\Li{[0,R]) < N\\g\\L^i^,uppr,)\\D'^f\\Lp(supp^), 

where N = N{p, m, K) > 0. 

Proof. In the case k = 1, the inequality (j7.7p is the classical Hardy's in- 
equality. For 1 < k < m, the inequality is known in one form or another. 
Here we giye a simple proof for completeness. By Taylor's formula, for any 
X G (0, R], we haye 

x-^D'^-^fix) = {V{k)Y^ / {l-rf-^{D'^f){rx)dr, 

Jo 

where T is the Gamma function. Thanks to Young's inequality, we get 

llx-'^I)— < imy' [\l-r)''-'\\{D^f){rx)U^^[o,R])dr 

Jo 

< {my' [\l-r)''-'r-'/^dr\\D^fU^^[oM 
Jo 

= r(i - i/p)(r(fc + 1 - i/p))"'||Z5"7llL,([o,ij]). 

The first assertion is proyed. The second assertion follows easily from the 
first one by using the Leibniz rule and Holder's inequality. □ 
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7.3. Estimates of mean oscillations. Now we prove the following esti- 
mate of mean oscillations. As in Section[6l we assume that all the lower-order 
coefficients of L are zero. 

Proposition 7.10. Let to e M, xq e R'l, Xq = {to,xo), r G (0,oo), k G 
[64, oo), A > 0, u e (2,oo), u' = 2u/{u-2), and = {fl . . . , fS^ G 
L2,ioc(M'^^)- Assume that nr < Rq and u G C;~(M!j+^) satisfies (fTIS]) and 

ut + {-irCu+Xu= ^"4 

\a\<m 

in Q'^^{Xq). Then under Assumption \2.'J[ (7^, there exists a function 
depending on Q+ := Q+,(Xo) such that N"^U <U^ < NU and 



where N = N{d, m, n, 5, u) > 0. 

The proof of the proposition is split into two cases. As in the proof of 
Lemma 15.21 we assume that coefficients aap and are infinitely differen- 
tiable. 

Case 1: the first coordinate of xq > Kr/16. In this case, we have 

Q+(Xo) = Qr{X^) C Q../16(^0) c \ 

Since k/16 > 4, ()7.8p follows immediately by applying Theorem 16. II with ^ 
in place of k. Note that Theorem 16. II is proved under the assumption that u 
vanishes outside Q-yRg. However, one can see that the proof of the theorem 
does not use this assumption in the case nr < Rq, more precisely, -^r < Rq. 

Case 2: < the first coordinate of xq < Kr/16. Denote yo = (0,Xg) and 
^0 = (to, Ho)- Without loss of generality, one may assume Yq = (0,0). Notice 
that in this case, 

QtiXo) C Q+ /g C Q+ C Q+ /2 C QtriXo)- (7.9) 

Denote R = Kr/2{< Rq). Because of Assumption 12.31 after an orthogonal 
transformation y = px centered at Yq = (0, 0), we may assume 

{{yi,y') : -fR < yi} n Br cnn Br c {{yi,y') : --fR < yijnBR, 

where O is the image of M^J. under the orthogonal transformation, and 

sup / \aap{t,y) -aai3{yi)\dxdt <-f\QR\. (7.10) 

Let Xq be the new coordinates of Xq after the orthogonal transformation. 
Then (|7.9|) becomes 

Cr(Xo) C Cr/4 C Cr/2 C Cr C CKriXo). (7.11) 
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Without any confusion, in the new coordinate system we stih denote the 
corresponding unknown function, the coefficients, and the data by u, aa/s, 
and /, respectively. Take a smooth function x defined on M such that 

Xivi) = for yi < jR, xivi) = 1 for yi > 2jR, 

\D^x\<N{nRy^ for A; = 1,2,..., m. 

Below we present a few lemmas, which should be read as parts of the 
proof of the second case. 

Lemma 7.11. Let u := xu. Then u along with all its derivatives vanishes 
on Qfi n {yi < 7-R} and satisfies in Q'J^ := Qr H {yi > ^R}, 

\a\ = \l3\=m 

+ xD^f^ + {-iTg + i-irh, (7.12) 

\a\<m 

where Cq is the differential operator with the coefficients Oap from ()7.10p . 
and 

g= J2 D''{a^pDPi{x-l)u)), 

\a\ = \l3\=m 

a| = |/3|=m 

Proof. This can be easily seen if one begins with multiplying the equation 
of M by X and then adding {—l)"^h to the both sides. □ 

Now let w be the unique if^iM x {y : yi > ^R}) solution of 

Wt + {-irCow + Xw=i-ir Y D''{^iaa(S-aai3)D^u^ 

\a\ = \l3\=m 

+ i-irg+ Y xD'^iM + (-l)'"^ (7.13) 

|a| <m 

in M X {y : yi > 7-R}, where (p is an infinitely differentiable function such 
that 

0<ip<l, ip = Ion Qji/2, ip = outside {-R^"", R^"") x Br, 

and 

g= J2 D"{aapipD^{{x-l)u)), 

\a\ = \l3\=m 

\a\ = \l3\=m 

Note that by the classical theory iv is infinitely differentiable in M x {y : 
yi > 7i?}. 
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Lemma 7.12. Let T G (-00,00] and fi' = {y G M"' : yi > 7/?}. Then for 
the function w in (17.13p . we have 

\a\<m \a\ = \(5\=m, 

\fS\=m |a|<m 

|/3|=m 

where N = N{d, m, n, 5) . 



+ ^ \\vD^'mL2{n'^n{{t,y)nB<yi<2'yR}), 



Proof. The first two terms on the right-hand side result from a direct appli- 
cation of Theorem 17.11 with ft' in place of M!^. For the third term, due to the 
presence of the factor X; one cannot directly apply Theorem 17.11 to get the 
desired estimate. However, we observe that as in the proof of Theorem 13.11 
after testing the equation by w and integrating by parts, the terms D^'u on 
the right-hand side of (jS.Sp are now replaced by D'^lxw). Notice that 

I-D'^xI < N{yi - jR)-^ for k = 1, 2, m. 

Then by Lemma 17.91 the L2 norm of D"{xw) is bounded above by the L2 
norm of D^'w, which implies that the same estimate as in Theorem 17. 1 1 still 
holds true even with the presence of x- For the last term in the above 
estimate, we argue in the same way and also use the property that x — 1 is 
supported on {y : yi < 2jR}. □ 

Recall that 

|o|<m |o|<m 

Lemma 7.13. For the function w in (j7.13p . we have 

m 

j;A^^(/Q.+ jZ?'=^P)y;<Ary/'^'([/^)y; + Ar(F2)y;, (7.14) 

fc=0 

where u and u' are from Proposition \7.10\ 

Proof. Note that ip{t,y) vanishes on {y : \y\ > R,yi > 'jR}. Thus by the 
estimate in Lemma 17.121 when T = 0, it follows that the left-hand side of 
()7.14p is less than a constant times 

\a\ = \l3\=m |/3|=m 



+ E >^^^H\mc'+ E {h.R<y.<^.R}\D'u\')^^ ■.= h+h+h+h. 

\a\<m |/3|=m 
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case we see that Ii is bounded by N^^/'^' {U'^)^^\ Observe that by Holder's 



It is clear that I3 is bounded by N{F . By using (|7.10p as in the interior 

case we see that Ii 
inequality we have 

{l{^R<yi<2iR}\D"''^\'^)c^ < {h'yR<yi<2iR})c^ (I^'^mDc^T 

<Njy^\\D"^un'J^. (7.15) 

Thus h is also bounded by N-/^^^' {U'')^^^ . 

To estimate I2, we notice that x ~ 1 = for t/i > 2'yR, D^~^x = in 
(-00, 7i?] n [27i?, +00), and 

\D^-Px\ < N{-iR)-\^\+\'^\ (7.16) 

in (7i?, 27/2), where = m and |/3| < m — 1. For any {s,jR,y') £ Qr, 
let iji = yi{s,y') be the largest number such that {s,yi,y') € dQ. (Indeed, 
in this case yi is uniquely determined as a function of y' .) Because yi € 
(—7/2,7/2), the inequality ()7.16p implies 

\D^-^x{s,y)\<N{y,-mr^^\^^^^- (7.17) 

We also notice that u, as a function of yi , vanishes along with its derivatives 
up to (m— l)-th order at (s, yi,y'). Thus, by LemmaEHtogether with (I7.17|) . 



r \D('{{x-l)u{s,yi,y'))\^dyi 



< \D^{{x-l)u{s,yi,y'))\^dyi<N \D^u{s,yi,y')\^ dyi, (7.18) 



yi -'yi 



where r = r{y') = min{27i2, R"^ — ly'p}- Integrating with respect to s 
and y' and using Holder's inequality as in ()7.15p . we estimate I2 by 

From this and the above estimates for Ii, i = 1, 3, 4, we conclude ()7.14p . □ 

Now we are ready to complete the proof of Proposition 17.101 We extend 
w to be zero in Cr \ Q]^ , so that w G H'^^Cr), and let w = + (1 — x)u. 
Since (1 — x)^^ vanishes for yi > 2jR, using the second inequality in (|7.18p 
and Holder's inequality as in (17.151) . we see that 

m 

A^^(|I)'=((1 - x)u)\%^; < 

k=0 

This combined with (|7.14p shows that 

{W%/^ < N^'/"' {Un'J; + NiF%/l (7.19) 
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Noting that Cr(Xo) C Cr and \CR\/\Cr{Xo)\ < N{d)K'^"'+'^ , we obtain from 

^^Xik) - {^'^"'(Uny: + {F%') . (7.20) 

Next, we define v = u — w in Cr. It is easily seen that v = in Cr \ Q]^ 
and V satisfies 

Vt + {-irCoV+Xv = Q 
in Qr/2 n {ui > 'yR} and vanishes along with its derivatives up to (m — l)-th 
order on Qr n {yi = jR}. Denote 

Vi=Cr{Xo)n{yi <-fR}, V2 = Cr{Xo)\Vi, V3 = QR/^r\{yi>-fR}. 

Because of (jT.lip . < NKj\Cr{XQ)\. Thanks to the fact that w is in- 
finitely differentiable in M x {y : yi > 7^2}, we see that v is infinitely 
differentiable in Qr n {yi > ^R}- Then applying Corollary 17.61 and Lemma 
17.71 with a scaling argument, we compute 

d^' - (^')c.(Xo)l)c.(Xo) + (I© - (0)c.(Xo)l)c,.(Xo) 
< Nr^/^{[V']ci/2^^^) + [e]ci/2(^,)) +iV^7ll^llL^(75.) 

which together with ([719]) and (f7:20|) yields (fT^Sj) . Indeed, we set U'^ = 
U' + e. Then 

(1^^ - (t^^)c,.(Xo)l)c.(Xo) ^ ^d^' - (^')c.(Xo)l)c.(Xo) 

+^(|0-(0)c,.(Xo)l)c.(Xo)+^Wc.(Xo)' 
where the last term is estimated by (j7.'2Up . As shown above, the first two 
terms on the right-hand side are estimated by N{k~^^'^ -\-K'y)(V'^)^^'^ ^, which 

is taken care of by (j7.19p and the fact that u = v + w in Cr. Finally, we 
transform the obtained inequality back to the original coordinates to get the 
inequality (j7.8p . This completes the proof of Proposition 17.101 

7.4. Proof of Theorem 12.41 We finish the proof of Theorem 12.41 in this 
subsection. As in Section [6l we assume all the lower-order coefficients of C 
are zero. 

Next in the measure space M."^^ endowed with the Borel ci-field and 
Lebesgue measure consider the filtration of dyadic parabolic cubes {Ci,l G 
Z}, where Z = {0, ±1, ±2, ...} and is the collection of cubes 

(io2-2-', {io + 1)2-2-'] ^ (.^2-', {h + 1)2-'] X ... X (irf2-', [i^ + 1)2-'], 

where iQ,ii, G Z, ii > 0. Notice that if X G C € C;, then for the 
smallest r > such that C C Qr{X) we have 

I I \g{Y)-g{Z)\dYdZ <N[d) I / \g{Y) - g[Z)\dY dZ. 
JcJc JQt{X)JQ+{X) 
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By using this, the following corollary is proved in the same manner as Corol- 
lary [Ol 

Corollary 7.14. Let 7 G (0, 1/4), A > 0, 1/ G (2, 00), i^' = 2v/{u - 2), and 

ssume that u G C^^(M.'^^) vanishes outside Q'^Ri^iZo) and 
satisfies (|7.3p as well as 

ut + {-ircu+xu= ^"/a 

\a\<m 



in M!^^ , where G L2joc(M'+~^) ■ Then under Assumption \2.3\ ("f), for each 
/ G Z, C G C;, and k > 64, there exists a function depending on C such 
that N-^U <U^ <NU and 

{U^-(U^')c\)^<N{F^)c, 
where N = N{d, 5, m, n, r) and 



Theorem 7.15. Let p G (2, 00), \ > 0, Zq e M^[+\ and G Lp(M^+^). 
There exist positive constants 7 G (0, 1/4) and N, depending only on d, 5, 

m, n, p, such that under Assumption \2.3\ ("f), for u G C'^^iM!^'^) vanishing 
outside Q-yniZo) and satisfying (j7.3p as well as 

ut + {-ircu+xu= ^"/a 

\a\<m 



in M'^^ , we have 



where N = N{d, 5, m, n, p) . 

Proof. It follows from the proof of Theorem 16.31 by using Corollary 17.141 and 
Theorem 13. 6i □ 



Proof of Theorem \2.4\ Due to the duality argument and Theorem 17. H it is 
enough to consider the case p > 2. In this case the theorem is proved using 
Theorem 1 7. 1 5 1 and the standard partition of unity argument. □ 

8. Systems on Reifenberg flat domains 

As an application, in this section we consider systems on a cylindrical 
domain (0, T) x Q. Here T G (0, 00) and is a Reifenberg flat domain in 
M'^, which is not necessarily to be bounded. Roughly speaking, the boundary 
of a Reifenberg flat domain is locally trapped in thin discs. 

We impose a similar regularity assumption on a^^ as Assumption 12.11 
Near the boundary, we require that in each small scale the direction in which 
the coefficients are only measurable coincides with the "normal" direction 
of a certain thin disc, which contains a portion of dQ. More precisely, we 
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assume the following, where the parameter 7 G (0, 1/20) will be determined 
later. 

Assumption 8.1 (7). There is a constant Rq € (0, 1] such that the following 
holds. 

(i) For any x G 17, t € M, and any r G (0, min{i?0) dist(x, dQ,)/2}] (so that 
Br{x) C ri), there is a spatial coordinate system depending on (f,x) and r 
such that in this new coordinate system, we have 

iai3{s,yi,y') - f aai3iT,yi,z')dz' dr dyds <-f. (8.1) 

(ii) For any x G 50, t G M, and any r G (0,i?o], there is a spatial 
coordinate system depending on [t, x) and r such that in this new coordinate 
system, we have (j8.ip and 

{{yi,y') : xi+7r < yi}r\Br{x) C Or(a:) C {{yi,y') : xi -7r < yi}nB,,(.x). 

We remark that the boundary of a Reifenberg flat domain may have a 
fractal structure. In particular, if the boundary is locally the graph of 
a Lipschitz continuous function with a small Lipschitz constant, then Vt is 
Reifenberg flat. Thus all domains are Reifenberg flat for any 7 > 0. 

The next theorem is the main result of this section. 

Theorem 8.2. Let T G (0, 00], ft be a domain in M'^, p G (l,oo), and 
fa = • • • 1 faY^ € Lp{flT), < "I- Then there exists a constant 7 = 
-y{d, n, m,p, 5) such that, under Assumption \8.1\ (7 ), the following hold true. 

(i) For any u G 'H^{il.T) satisfying 

ut + {-irCu+Xu= D'^f^ (8.2) 

\a\<m 

in r^T, we have 

^ X^-^\\D^u\\L^^nr)<N E ^^ll/allL,(n,), (8-3) 

|o|<m |ci|<m 

provided that A > Aq, where N and Aq > depend only on d, n, m, p, 5, K , 
and Rq . 

(ii) For any A > Aq, there exists a unique u G %'^{VLt) satisfying ()8.2p . 
(Hi) Let A = and suppose T G (0,00). There exists a unique solution 
u G ni\{Q,T) X fl) of ([O]) with the initial condition ii(0, •) = d in Q.. 
Moreover, u satisfies 

ll^lk™((0,T)xf7) < ^ E ll/allip({0,T)xn), 
|ci| <m 

where N depends only on d, n, m, p, 6, K , Rq, and T. 

A similar result for second-order scalar elliptic equations with symmetric 
coefficient matrices in bounded Reifenberg flat domains was recently stud- 
ied by Byun and Wang [6]. The proofs there are based on an approach 
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developed by the same authors in which in turn uses an idea of approx- 
imations originally due to Caffarelli; see [8]. This approach is quite suitable 
for studying equations on domains with rough boundaries. Here we adapt 
it to investigate higher-order systems. Our proofs, however, are in several 
aspects different from those in [6]. In particular, we do not use the reverse 
Holder's inequality, the iteration argument or the maximum principle, the 
latter of which is not available in our case. Compared to [6], we consider 
more general operators and domains by allowing lower-order terms and un- 
bounded domains. 

The crucial ingredients of the proofs below are the interior and the bound- 
ary estimates established in Sections U] and [71 By a scaling, we may assume 
i?o = 1 in the sequel. Recall the definitions of U, V, W, and F in Sections 
land El 

Lemma 8.3. Let R G (0,1], A G (0,oo), z/ G (2,oo), u' = 2u/{v - 2), 
fa = ifL • • • > /a)**^ ^ L2jocO^ X i}) , \a\ < m. Assume that = for any 
a, (3 satisfying \a\ + < 2m and that u G C^(M x 0) satisfies ()8.2p . Then 
the following hold true. 

(i) Suppose G 0, dist{0,dfl) > R, and Assumption \8.1\ ('y) (i) holds at 
the origin. Then u admits a decomposition u = v + w in Qn, and w and v 
satisfy 

{W')'£ < Njy-^'iUn'd'^ + N{F')'^l (8.4) 

and 

WnL^iQn;.) < ^^'""'{UnZ + NiF^)Tn + (8.5) 
where N = N{d, n, m, 5,v) > Q is a constant. 

(a) Suppose G dVL and Assumption \8.1\ (ii) holds at the origin. Then 
u admits a decomposition u = v + w in Cr, and w and v satisfy 

and 

WVh^^Cn,,) < N^'^^'iUn'J: + N{F^)'j;^ + iV([/2)J/^ (8.7) 
where N = N{d, n, m, 6,u) > is a constant. 

Proof. The proof is similar to that of Proposition 17.101 with some modifica- 
tions. Without loss of generality, we may assume Assumption 18.11 holds in 
the original (i, x)-coordinates. First we prove Assertion (i). Define 

«Q/3(yi) = / aai3{T,yi,z')dz' dr, 

and Cq be the operator with in place of Oq,^ in C. Then in Q^, u satisfies 

ut + {-irCQU+\u 

l^^l^l/^l^"^ |Q|<m 
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Now let w be the unique 

^m(]gd+i) solution of 
Wt + {-l)"'CQW+\w 

= |a|<7?i 

in M'^^^, where ip is an infinitely differentiable function such that 

0<ip<l, 99 = 1 on Qr/2, ip = outside (-i?^", R^"") x Sj^. 

The existence of such solution is due to Theorem 13. II By the same theorem, 
we have 



\a\<m 

+N ^ ^^ip{aai3 - aai3)D'^u 

\a\ = \l3\=m 

where Mq = (— oo, 0) x M'^. This together with Holder's inequality gives 
(j8.4p . Next, it is easily seen that v := u — w satisfies in Qji/2 

Vt + {-irCov + Xv = 0. (8.9) 

With a scaling, we apply Corollary 14.71 to (j8.9p to get 

The estimate above together with (j8.4p leads to (j8.5p . 

Next, we prove Assertion (ii). Define Oq,^ as above and notice that u 
satisfies (|8.8p in Cr. To apply Corollarv 17.81 we need to locally cutoff u so 
that the new function vanishes for yi < 'jR. To this end, take the function 
X from the proof of Proposition 17.101 Then ii := along with all its 
derivatives vanishes on yi < 'yR and satisfies ()7.12p in Qj^ := Qj?, n {yi > 
jR}. Let w be the unique HfiM. x {y : yi > jR}) solution of (|7.13p in 
M X {y : yi > -jR}. Following the argument in the proof of Proposition 17. lOl 
with obvious modifications, we get 

m 

^A^^(/Q..|Z?^-H^C<A^y/^'(f/^)y; + A^(F^)y,^ (8.10) 

k=0 

We extend w to be zero in C/j \ Q]^ , so that w G 'H™(Cr), and let w = 
w+ (1 — x)u. By Lemma 17.91 and Holder's inequality, we get (|8.6p from 
(I8T0]) . 

Next, we define v = u — w in Cr. It is easily seen that v = in Cj^ \ Q^j^ 
and V satisfies 

vt + {-ir Cqv + \v = Q 

in Qfi/2 n {yi > ^R} and vanishes along with its derivatives up to {m — l)-th 
order on Qji n {yi = jR}- Applying Corollarv 17.81 we get 

ll^llLoo(Cfl/4) = ll^llLoo(Qfl/4n{?/i>7R}) -^(^^)ci^/2' 
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which together with ()8.6p gives ()8.7p . This completes the proof of the lemma. 

□ 

For a fmiction / on a set T? C M"^"^^, we define its maximal function A4f 
by Mf = A4(/p/). For any s > 0, we introduce two level sets 

A{s) = {{t,x) £Rxn:U>s}, 

B{s) = {(i,x) € M X : 7-i/'^'(A^(f2))^/^ + {MiU"))^^" > s}. 
With Lemma 18.31 in hand, we get the following corollary. 

Corollary 8.4. Under the assumptions of Lemma \8.3\ suppose & Cl and 
Assumption \8.1\ ('y) holds. Let s S (0, oo) be a constant. Then there exists 
a constant n € (l,oo), depending only on d, n, m, 5, and v, such that the 
following holds. If 

|C«/32n^M| >7'/"'|Ck/32|, (8.11) 

then we have C B{s). 

Proof. By dividing u and / by s, we may assume s = 1. We prove by 
contradiction. Suppose at a point (t, x) S Cj:j/32, we have 

^-^Z"' {M{F^){t,x)f'^ + {M{U^){t,x)f''' < 1. (8.12) 

Let us consider two cases. 

Case 1: dist{0,dn) > R/8. Notice that 

{t,x) e Cr/32 = Qr/32 C Qr/s cRxn. 

Due to Lemma [831 (i), we can write u= w+ v in Qj^/s ^i^d, by ()8.12p . 

where A'^i and constants Ni below depend only on d, n, m, 6, and u. By 
()8.13p . the triangle inequality and Chebyshev's inequality, we get 

|Ci?/32n^(K)| = \{{t,x) G : [/ > k}| 

< \{it,x) e Cr/32 :W>K-N,}\<iK- Ni)-^N^j^/'^'\Qn/,\, 

which contradicts with ()8.1ip if we choose k sufficiently large. 

Case 2: dist{0,d^l) < R/8. We take y £ dQ such that \y\ = dist(0,9O). 
Notice that in this case we have 

{t,x) G C C^/4(0,y) C Cjf(0,y). 

Due to Lemma [8^ (ii), we can write u = w + v in Cr{0, y) and, by (I8.12|) . 

II^IIl^(c«/,(o,.)) < N,, (W^^)i2o,,) < N,^'/'^'. (8.14) 
By (j8.14p . the triangle inequality and Chebyshev's inequality, we get 

\CR/32nA{K)\ = \{{t,x) G Cr/32 : [/ > k}| 

< \{{t,x) G C^/32 ■.W>k-N2}\<{k- N2r^Nh^/'''\CR\, 
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which contradicts with (18.111) if we choose k sufficiently large. □ 

Theorem 8.5. Let p € (2,oo), X > 0, Xq e W^^'^ and G Lp{R x O). 
Suppose that aap = for any a, 13 satisfying \a\ + \(3\ < 2m, and u G 
C^(]R X Q) vanishes outside Q^{Xq) and satisfies (j8.2p m R x $7. There 
exist positive constants 7 € (0, 1/20) and N , depending only on d,6,m,n,p, 
such that, under Assumption \8.1\ (^) we have 

\\U\\Lp(Rxn) < ^l|-^llLp(I8xn)> (8.15) 

where N = N(d, 5, m, n, p) . 

Proof. We fix = p/2 + 1 and let v' = 2v/{u — 2). Let k be the constant in 
Corollary 18. 4[ For any s > 0, by Chebyshev's inequality, 

\A{ks)\ < i^syiml^^^^^y (8.16) 

From (j8.16p , Corollary 18. 4^ and a result from measure theory on the "crawl- 
ing of ink spots" which can be found in [36] or [26^ Sect. 2], we have the 
following upper bound of the distribution of U (cf. [U Sect. 4]). For any 
7 G (0,1/20] and tis > 7"'/"' llf^llL^CExn), 

\A{ks)\<N^^^/'''\B{s)\. (8.17) 

For < Ks < T^"*^^"^' ||?7||£,2(]Rxn), we simply bound the distribution function 
by (I8.16|) . Now we use the elementary identity: 

> / \{{t,x)eV:\f{t,x)\>s}\sP-Us, 
Jo 

to deduce from (KTEh and (fHTT]) that 

By Holder's inequality, 

ll^lliaCKxn) = \\U\\L2{Q-.f{Xo)nRxn) < ^||f^||Lp(MxQ)7^'^'^^"*^^^^^~"^^^^ • (8.18) 
Since 2 < u < p, hy the Hardy-Littlewood maximal function theorem and 
(f8J8]) . we obtain 

iir/ii^ < Arfi'v(2-p)/i-'ii 77||P + ATfi'Y^/'^'iir/iiP 

llLp(Mxf7) - ^^67 11-^ llLp(Rxn) ^ ^^67 II ^ llLp(Kxn) • 

To get the estimate (j8.15p . it suffices to take 7 = 7(0?, n, m, (5, p) G (0, 1/20] 
sufficiently small such that Na-f'^/"' < 1/2. □ 

Now we are ready to give the proof of Theorem 18. 2i 



p 



Proof of Theorem \8.S\ , We derive Assertion (i) from Theorem 18.51 as in the 
proof of Theorem 12. 2i Assertion (ii) follows from Assertion (i) , the method 
of continuity and an approximation argument. To prove Assertion (iii), let 
u = e*(^o+^)M. Clearly u satisfies with Aq + 1 and := e-*(^o+^)^ in 
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place of A (= 0) and / respectively. Therefore, to show the existence and 
uniqueness of u is equivalent to show those of the solution ii to the equation 

Vt + {-ir£v+{Xo + l)v= ^"/a (8-19) 

with the zero initial condition. We extend / to be zero for t < 0. Using 
Assertion (ii), we find a unique T-1'^{Qt) solution v of (IS.lOh in Ut- Set 
u := V. By the uniqueness, v = for t < 0, which implies that ii satisfies 
the same equation with the zero initial condition at t = 0. This gives the 
existence. For the uniqueness, let ii G ?i™((0, T)r2) be a solution to (j8.19p 
with the zero initial condition and zero right-hand side. We extend ii to be 
zero for t < and denote it by v. It is easily seen that v satisfies (jS.lOp 
in ^It with zero right-hand side. By Assertion (ii), we have t; = in 0,t- 
Finally, the estimate of u is deduced from the estimate of u, which in turn 
follows from Assertion (ii). The theorem is proved. □ 

We finish the paper by considering the corresponding Dirichlet boundary 
value problem for elliptic equations in the case that all the involved functions 
are independent of the time variable. 

Theorem 8.6. Let he a domain in , p G (1, oo), andf^ = (/^, . . . , /„ )*'' G 
Lp(Q), \a\ < m. Then there exists a constant 7 = 'j{d,n,m,p,6) such that, 
under Assumption \8. 1\ ('jj, the following hold true. 

(i) For any u G W^{Vl) satisfying 

Cu+{-l)"'\u= D'^f^ in n, (8.20) 

\a\<m 

we have 

Y X'-^\\D^u\\L^^n)<N Y A^||/J|L,(n), 

|a|<m |o|<m 

provided that A > Aq, where N and Aq > depend only on d, n, m, p, 5, K , 
and Rq. 

(ii) For any A > Aq, there exists a unique u G Wp'-{^) satisfying (|8.20p . 
(Hi) Let A = 0. Assume that < 00 and that equation ()8.20p admits a 

estimate, i.e. for any u G W2^{Q) satisfying (I8.20p with G L2{i^), 
we have 

\\u\\w^in) Wfah^in), 

\a\<m 

where N depends only on d, n, m, 5, K, Rq, and Then there exists a 
unique solution u G W'^{^) of (|8.20p and u satisfies 

\\u\\w^^^)<N Y ll/allL,(n), (8.21) 

|a|<m 

where N depends only on d, n, m, p, 6, K, Rq, and 
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Proof. The proofs of the first two assertions in Theorem 18.61 are completely 
analogous to those of Theorem 18.21 and are thus omitted. 

For the last assertion, as before it is sufficient to show the a priori estimate 
(j8.2ip for u € Wj7'(^)- Again, by a duality argument we may focus on the 
case p > 2. From Assertion (i), we have 

M\w,^{n) <N Wfah.in) + N-r\\u\\L^^n), (8.22) 

\a\<m 

Therefore, we only need to estimate Takepi € (^,00) such that 1 — 

d/p > —d/pi. By Holder's inequality. Young's inequality and the Poincare- 
Sobolev inequality, we get for any e > 0, 

||tt||Lp(n) < A^(e)l|ii||L2(n) +e||it||Lp,(f7) < N{e)\\u\\L^(^^) + Niie\\Du\\L^(^^y 

Choosing e = l/{2NjN^) and using ([522]), we obtain The theorem 

is proved. □ 

Remark 8.7. As an immediate corollary of Theorem 18.61 we obtain the 
W^{9) solvability of (f8:20|) with A = if, in addition, < 00 and the 
lower-order coefficients of C are all zero. Finally, for second-order scalar 
elliptic equations in the form 

Dj{aijDiu) + Di{aiu) + hiDiU + cu = div (7 -|- /, 

under the same assumptions and the additional condition that DjOj -|- c < 
in VL in the weak sense, we also get the Wj^(fl) solvability by using the 
classical W2 estimate of the Dirichlet problem on a domain with a bounded 
measure; see |20l §8] or \\.4\ Sect. 7]. This result generalizes Theorem 2.7 
in |14] . in which bounded Lipschitz domains are considered. It also extends 
the main result of [6] to equations with lower-order terms. 
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